Pairing with Unconventional Symmetry around BCS-BEC Crossover: Fermionic 
Atoms in 2D Optical Lattices to Correlated Electron Systems 



Dusan Volcko and Khandker F. Quader 
Department of Physics, Kent State University, Kent, OH J^J^2^2 
(Dated: February 2, 2008) 

We study superfluid properties of fermions on a 2D lattice using a finite-range pairing interaction 
derivable from an extended Hubbard model. We obtain signatures of unconventional pair-symmetry 
states, d^2_y2 and extended-s (s*), in the BCS-BEC crossover region. The fermion momentum 
distribution function, v1, the ratio of the Bogoliubov coefficients, Vk/u]^, and the Fourier transform 
of are among the properties that are strikingly different for d- and s* symmetries in the crossover 
region. Fermionic atoms in 2D optical lattices may provide a way to observe these signatures. We 
discuss possible experimental ramifications of our results. 
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Attainment of boson and fermion condensates [l| of 
ultracold neutral atoms has presented an unprecedented 
opportunity to study properties of quantum many- 
particle systems. Fermionic atoms in optical lattices 0,0] 
constitute yet another intriguing set of systems. While 
these are by themselves interesting to study, they may 
also provide a way to gain useful insight into proper- 
ties of correlated electrons in solids. Jaksch et al [i\ 
suggested that atoms in optical lattices, confined to the 
lowest Bloch band, can be represented by the Hubbard 
model with hopping kinetic energy t between neighboring 
sites, and on-site interaction U . Hubbard model calcula- 
tions P, @, 01 predict that attractive-U Hubbard model 
give rise to s-wave superconductivity, while the repulsive- 
U model results in an antiferromagnetic or a d-wave 
superconducting phase depending on filling (number of 
fermions per lattice site). Owing to the continuous tun- 
ability of model parameters such as, density, hopping or 
interactions, optical lattices can serve as testing grounds 
for such models. This has led, for example, to the sug- 
gestion 0] that the underlying physics of the high Tc 
superconductors may be understood by studying these 
systems. Recent work have pointed out possible 

role of additional Bloch bands and multi-band couplings 
in optical lattices. In solids this would correspond to 
having multiple orbitals and near-neighbor interactions. 
Duan Q has shown that on different sides of a broad 
Feshbach resonance, the effective Hamiltonian can be re- 
duced to a t-J model, familiar in correlated electron sys- 
tems, wherein it has been suggested that t-J model 
can give rise to d-wave pairing. 

Fermionic atoms subjected to positive and negative de- 
tuning using Feshbach resonance technique provide real- 
izations of BEC-BCS crossover behavior. It has been 
recently suggested [ll| that it should also be possible 
to study superfluid properties of fermions in optical lat- 
tices around BEC-BCS crossover regime. Starting with 
the seminal work of Eagles and Leggett the 
BEC-BCS crossover problem received considerable the- 
oretical attention d, H, [H, [H, [13, [11, [11 due to 



the possibihty that high Tc superconductors, possess- 
ing short coherence lengths, could fall in the BEC-BCS 
crossover region. Several authors employed continuum 
models [ij, [l^, [H, [l^ j focussing mostly on conven- 
tional s-wave pair symmetry. Lattice models with on- 
site or nearest- neighbor attractions have also been con- 
sidered [i,[Il[I3,[ll[ll- More recent theory work ^ 
are in the context of cold fermions. 

Motivated by these issues, in this paper, we study su- 
perfluid properties of fermions in a 2D square lattice in 
the BEC-BCS crossover regime using a finite-range pair- 
ing interaction, obtainable from a multi-band extended 
Hubbard model. As representative cases of unconven- 
tional pair symmetry, we consider two even-parity repre- 
sentations of the cubic group, namely the £ = 2 dx2_y2- 
wave, and the £ = extended s-wave (s*). There has 
been work [isl [Til [2]| employing similar pairing inter- 
action; however these have focussed on different systems 
and issues. We present several new results, including spe- 
cific signatures of superfluid states with unconventional 
pairing gap symmetry as one goes between the BEC and 
BCS regimes. This could provide a way to distinguish be- 
tween different gap symmetry states in systems that al- 
low for tuning into the BEC-BCS crossover regime, such 
as fermionic atoms in 2D optical lattices, and possibly 
high Tc cuprates. One of our key results is the remark- 
able behavior of the fermion distribution function, w^, 
(related to momentum distribution, n^): For the d-wave 
gap function, v'^ changes abruptly from having a peak at 
the Brilloiun zone (BZ) center (0,0) to a vanishing cen- 
tral peak accompanied by a redistribution of the weight 
around other parts of the BZ ((0, ±7r),(±7r, 0)) as the sys- 
tem crosses from the weak-coupling BCS to the strong- 
coupling BEC regime. By contrast, changes smoothly 
in the s*-wave case. Similar signatures are also found in 
the ratio of Bogoliubov coefficients Vk/uk, related to the 
phase of the superfluid wavefunction. The Fourier trans- 
form of I'l in real space exhibits a "checkerboard" type 
pattern that could have consequences for experiments. 

The extended Hubbard model for two equal species 
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population system on a 2D square lattice is given by: 
H = (-^vCj<,+H.c.) + U^ni<,ni_<, 

<ij><y i 
~ ^ X! "io-njcr' - Mo^"-j, (1) 

where t is the kinetic energy hopping, /io the unrenormal- 
ized chemical potential, U the on-site repulsion and V the 
nearest-neighbor attraction. In the case of cold fermions 
on a lattice, V would be related to inter-band coupling. 
a is the "pseudo-spin" index, that could refer to equally 
populated hyperfine states in the case of optical lattices. 
At the mean-field level, the Hartree self-energy terms 
renormalize such that /i = /.to + l-J-uif) + P-vif) where 
fJ-uif) and fJ-vif) are filling-dependent corrections to fi. 
We work with the renormalized fi so as to properly deal 
with weak and strong couplings, and take fij^if) = Jif, 
where = U, -V. The filhng / = N/2M, with N the 
number of particles, M the number of lattice sites, and 
the pseudo-spin degeneracy factor 2. On Fourier trans- 
forming and retaining interactions between particles with 
equal and opposite momentum, as in BCS theory, the re- 
duced pairing Hamiltonian assumes the form: 

Hpair ^Y^{<^k - fJ')ctck + ^ Vfcfc' C^, C+ j,, C-fcCfe (2) 
k kk' 

where in the tight-binding approximation, — 
-2t{coskx+cos ky); Vuk' = Vo{cos{k^-k',j.)+cos{ky~k'y)), 
which is non- separable. Using the standard BCS vari- 
ational ansatz, \^bcs >= HkC^k + i'kcjjcLk)|0 >, we 
obtain the T = gap equations for the gap functions 
Aj!'* = Ao(/)(cosfca; ibcosfcy) with dx2_y2{-) and s*{+) 
symmetries, 

1 1 cosfc2;(cosfca: ± cosA:y) 

where = ((efc - /i)^ + A2(cos /c^ ± cos fcj,)2)i/2. The 

Bogoliubov coefficients are given by, 

l-kp; KP = ^(i±^)- (4) 

The ratio v^uy, = -{E"^''' - {ck - /.t))/Af'"*. Follow- 
ing Leggett[13], we readjust /x for strong attractions by 
supplementing the T=0 gap equation with the number 
equation: 

BZ 

This determines the self-consistently readjusted /i, which 
is no longer fixed at the Fermi level, and makes the 
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FIG. 1: Chemical potential ^ vs. coupling V at different 
fillings / for the d-wave case. BEC pairs appear where iJi.{V) 
crosses the /i/2t = —2 line. The inset shows ^{V) for the s- 
(dash-short dashed line), s* (dashed line), and d-wave (solid 
line) at / = 0.2. 



gap equation applicable over the entire range of filling, 
thereby the BCS and BEC regimes. To allow for strong 
scattering, the sums arc performed over the entire BZ. 
The natural momentum cut-off afforded by the lattice 
avoids any possible ultraviolet divergences. 

Remarkable differences in features stem in an essential 
way from differences in gap symmetry. The (1^2 _y2 gap 
AJ? vanishes along the lines ±kx = ±fc.y in the 2D BZ, 
i.e. at four points on the Fermi surface(fs), the location 
of which depends upon filling. The s* gap A^ coincides 
with the tight-binding fs at exact 1/2- filling, and is node- 
less otherwise. Here, < 0, with /i = — 4i at the bottom 
of the band. Owing to particle-hole symmetry, it is suffi- 
cient to consider < / < 1/2. Upon examination of the 
gap functions and Eqs. (3-5), the following distinctions 
become apparent: 

(a) For very low fillings (/ — > 0, /i — 4t), a thresh- 
old coupling is required for pairing in the d-wave case, 
while in the s* case A'' — > as V ^ due to a weak 
singularity at fi = —At. On the other hand, at 1/2-filling, 
due to a weak singularity at /i = in the d-wave case , 
A'* — + as y ^ 0. In the s* case such a singularity is not 
present and as A* 0, V/4t tt^/S, i.e. a minimum 
coupling is needed for pairing. In contrast with A^ (V), 
A'^{V) changes slope at ^ — —At, and hence not smooth 
everywhere (though continuous). 

(b) For small fc, we have the following limiting be- 
havior: (i) efc < /i(= -At); \uk\ l,|wfc| 0; 
this is the strong- coupling BEC limit. Here the ratio 
Vk/uk Afc/2|/x| — > {kl - fc^)/2|/i|, i.e. analytic, (u) 
tk > l^i— —At); \uk\ 0,\vk\ 1; this is the weak- 



3 





1 


' 1 ' 


2.5 




d-wave 


2 


s -wave 




1.5 




/ - y / 


1.0 








,/^^^ 


/// / 


0.5 


/// / 







/ / / / 




0.5 1.0 1.5 2.0 2.5 3.0 

/ / 


/// / 

' / / 
/ / 




■ ■ / / 
/ ' / 

-' , / 


/ / 
/ M=-4-t 






/ /i=-3t 




' 1 


' {L=-Z\. 




/ / / / 


- -11=0 
















1 


< 1 , 



12 3 

V/4t 



FIG. 2: d-wave gap functions A/2t vs. nearest-neighbor cou- 
pling V/4:t for different chemical potential jj,. Inset : Results 
for the s* case, fi = —it demarcates BEG and BCS regimes. 

coupling BCS limit. Here Vk/uk — > l/{kx — ky), i.e. non- 
analytic, (iii) efc = /i(= — 4t); \uk\ ^ 0, \vk\ ^ 0, when 
Ek 0. Then Vk/uk {kx - ky)/{kx + ky), i.e. inter- 
mediate between (i) and (ii). It may be noted that for 
d-wave, the quasiparticle excitations in the BCS limit (ii) 
are "gapless" for some values of k, while in the BEC limit 
(i), Ek 0, even for gaps with nodes [2^ . 

Self-consistent numerical solutions of Eqs.(3-5) bear 
out the above features in detail, and also reveal a number 
of other features. We scale fi, V, A by hopping parame- 
ter, t. At a given fiUing /, both Af and A| increase with 
increasing V. While for d-wave it is easier to pair elec- 
trons at higher fillings, this is not necessarily the case for 
s*-wave for the weaker couplings V/At < 1.5 and small 
gaps A'' /2t < 0.5. 

In Fig.l we show ^{V) for different fillings f. At a fixed 
/, in both the d- and s*-wave cases, decreases with 
increasing coupling V, changing less rapidly for progres- 
sively larger /. However in the s* case, Ijl{V) exhibits a 
small "bump" for weaker couplings V/At < 1.5. The drop 
in fj, with increasing attraction is significantly more rapid 
in the uniform s-wave case; see inset in Fig. 1. Crossover 
to the BEC regime here is signalled by fJ-{V) going below 
the bottom of the band, i.e. crossing the fi — —At line. 
As Fig. 1 shows, for the d-wave case, this develops at 
both low and high fillings at some minimum value Vb/At 
of the coupling. It is interesting to note that as / — * 0, 
Vf,/4< 1.8. At exactly 1/2-filling this coupling tends 
to infinitely large values. For couplings F > 14, the sys- 
tem is conducive to BEC pairing; for V <Vb, the system 
exhibits BCS-like features. 

Fig. 2 shows the behavior of the d-wave gaps as a 
function of coupling V for different values of the chemi- 
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FIG. 3: (a), (b): 3D plots of d-wave electron distribution 
functions v'f. vs. kx — ky at filling / = 0.1, showing abrupt 
"jump" in vl. In BCS regime (a), ^l = -3t, A'' = 5.2t, 

V = 18. 7t, and in BE regime (b), ^l = -6t, A"^ = 0.6t, 

V = 5.2t. (c), (d): 3D plots of d-wave Vk/uk vs kx — ky 
for the same parameters as in ((a),(b)) respectively. In BCS 
regime (c) it can be seen to be non-analytic; in BEC regime 
(d) it is analytic, (e), (f): The same as in (a), (b), but for 
s*-wave; the behavior is smooth. 



cal potential ii. The /i = —At curve represents the locus 
of Vb/At for different fillings (see Fig. 1), and demarcates 
BEC and BCS -pair regimes. To the left is the fj, > —At 
region wherein finite gaps of the BCS or intermediate 
BCS-BEC types exist. On a given constant-// curve it 
may not be possible to have solutions for any arbitrary 
filling, but only those that satisfy Eqs. (3) and (4) self- 
consistently. The inset in Fig. 2 shows the corresponding 
A* (y) curves for the s* case. There are interesting dif- 
ferences with the d-wave results in that the boundary 
(/i = —At) separating BEC/BCS regimes is not as clear- 
cut for the weaker couplings V/At < 1.5 and the smaller 
gaps A/2t < 0.5, however the /i < —At region lies to the 
right of the fi — —At curve as in the d-wave case. 

Differences in the gap symmetry manifest in a strik- 
ing manner in the momentum distribution function, v^, 
and the ratio Vk/uk- For d-wave, for a given filling, in 
the weak-coupling BCS regime {V < Vb{f), fi > —At), 
vl exhibits a peak centered around the zone center (0,0), 
that becomes progressively narrower with decreasing fill- 
ing. Then at the crossover point at Vb{f) (/U = —At), vi 
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FIG. 4: (a) Fourier transform pv{x,y) of a typical d-wave 
electron distribution function, v^. Here, filling / — 0.01, fi — 
—4:.2t (strong-coupling regime), gap A = .76t. (b) Projection 
of (a) to show contrast ratio of pv{x,y). 



abruptly goes to zero around (0,0), and shows a drastic re- 
distribution in a different region of the BZ, namely, along 
(0, ±7r), and (±7r, 0). The abruptness is evident from the 
"jump" in as the chemical potential goes from just 
above the bottom of the band (/i > — 4t) to just below 
(/Lt < — 4t), i.e. from BCS to BEC regime. A representa- 
tive case is shown in Figs. 3a, 3b. In marked contrast, in 
the the s*-wave case (Figs 3e,3f), the zone center peak in 
decreases smoothly as one goes from the BCS regime 
to the BEC regime; only a slight redistribution occurs at 
(±7r, ±7r). We find this behavior to be replicated at all 
fillings f. As observed above in the limiting cases, the nu- 
merical calculations show (Fig 3c, 3d) that for d-wave, in 
the weak-coupling BCS regime, Vk/uk is non-analytic at 
±kx = ±fcj/; in the strong- coupling BEC regime, Vk/uk is 
analytic, vanishing along the zone diagonals and peaking 
about (±7r, 0), (0, ±7r). In the s* case (not shown), Vk/uk 
is analytic in both regimes. Similar behavior in Uk has 
also been reported in other work 13 ■ 

Our findings suggest that experiments, that may be 
able to directly or indirectly probe or combinations of 
Uk and Vk, could reveal novel aspects of the paired states. 
For example, it may be possible to decipher the OP sym- 
metry (e.g. d- or s*- wave) by measuring v'^ as a function 
of filling (especially at low-fillings), and/or for different 
interaction strengths, both of which can be controlled in 
optical lattices. At the BCS-BEC crossover, we expect 
the behavior to be quite different depending on whether 
the OP is d- or s* wave. Also, in the case of d-wave pairs, 



quantities sensitive to or to (ufc, Vk) should be very dif- 
ferent depending on whether the paired state is BEC or 
BCS like. A possible probe may be ARPES. Information 
may also be obtained from experiments that sample the 
quasiparticle energy E'fc = (A^.-|-e^)-^/^ (related to m^, Wfc), 
the quasiparticle density of states, or coherence factors, 
UkVk + u'i.v i,. A ngle-dependent or transverse ultrasound 
attenuation [23|, or quasiparticle tunneling at low fillings 
are possible experiments. 

The Fourier transform of v^{kx,ky), namely, pv{x,y) 
may provide yet another interesting way to test our 
results. In the d-wave case, in marked contrast with 
its behavior in the BCS regime, py{x,y) is oscilla- 
tory in the BEC regime, and exhibits an inhomogenous 
"checkerboard-type" pattern as shown in Fig 4(a,b). For 
the chosen parameters of Fig. 3, the contrast ratio of the 
lowest density to the peak is roughly 50%, being most 
sensitive to the location of pi{V). The length scale is of 
the order of fractions of lattice spacing. py{x,y) is fairly 
uniform in the s* case in both regimes. Highly sensitive 
STM may be able to pick up such distinctions [IS] • 

Much of the phenomena we have discussed are away 
from exact 1/2-filling, and at relatively strong coupling, 
where possible effects of spin density wave (SDW) and 
charge density wave (CDW) instabilities are expected to 
be suppressed. Addition of a next-near-neighbor hop- 
ping would also stabilize the paired state, as well as 
lower the minimum near-neighbor interaction necessary 
for a bound-state; we have checked this^M]- We have 
not explored here the issues of collective modes or phase 
separation ^25). It may be interesting to extend this work 
to, for example, finite-T, or to explore whether the inho- 
mogenous density that we find bear relationship to the 
range/strength of the interaction, or to possible phase 
separation. 
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